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ABSTRACT 

Aims. The Herschel survey of the Galactic Plane (Hi-GAL) provides a unique opportunity to study star formation over large areas of the sky and 
different environments in the Milky Way. We use the best studied Hi-GAL fields to date, two 2° • 2° tiles centered on it, b) = (30°, 0°) and (I, b) 
= (59°, 0°), to study the star formation activity in these regions of the sky using a large sample of well selected young stellar objects (YSOs). 
Methods. We use the Science Demonstration Phase Hi-GAL fields, where a tremendous effort has been made to identify the newly formed 
stars and to derive as accurately as possible their properties such as distance, bolometric luminosity, envelope mass and stage of evolution. We 
estimate the star formation rate (SFR) for these fields using the number of candidate YSOs and their average time scale to reach the Zero Age 
Main Sequence, and compare it with the rate estimated using their integrated luminosity at 70 /jm combined with an extragalactic star formation 
indicator. 

Results. We measure a SFR of (9.5 ± 4.3) • 10~ 4 M G /yr and (1.6 ± 0.7) • 10~ 4 M G /yr with the source counting method, in t = 30° and t = 59°, 
respectively. Results with the 70 /jm estimator are (2.4 ± 0.4) ■ 10~ 4 Ms/yr and (2.6 ± 1.1) • 10~ 6 M Q /yr. Since the 70 /jm indicator is de- 
rived from averaging extragalactic star forming complexes, we perform an extrapolation of these values to the whole Milky Way and obtain 
SFR Mlv = (0.71 ± 0.13) M Q /yr from 1=30° and SFR MH ,= (0.10 ± 0.04) M Q /yr from t = 59°. The estimates in t = 30° are in agreement with the 
most recent results on the Galactic star formation activity. 

Conclusions. The source counting method gives results valid only for the particular region under consideration. On the contrary, the IR indi- 
cator, by construction, gives results which can be extrapolated to the whole Galaxy. In particular, when the extragalactic indicator is applied to 
the £ = 30° field, it provides a SFR consistent with previous estimates, indicating that the characteristics of this field are likely close to those of 
the star-formation dominated galaxies used for its derivation. Since the sky coverage is limited, this analysis will improve when the full Hi-GAL 
survey will be available. It will cover the whole Galactic Plane, sampling almost the totality of Galactic star forming complexes. By means of 
the candidate YSOs counting method it will be then possible to calibrate a SFR Galactic indicator and to test the validity of the extragalactic 
estimators. 

Key words. Stars: formation, Galaxy: stellar content, Surveys 



1. Introduction 

The current estimates of the Star Formation Rate (SFR) of the 
Milky Way (MW) are uncertain mainly because we lack a gen- 
eral knowledge of the structure of our Galaxy. The most recent 
data suggest it is a two-armed, barred spiral with several sec- 
ondary arms bu t the actual number and position of the arm s 
is still unclear [Dameetal., 2001; Dame & Thaddeus, 1201 in . 
Moreover, the location of the Solar System in the Galactic Plane 
makes the definition of the overall Galactic structure, and conse- 
quently the heliocentric distance determination of star forming 
regions, even more difficult. Observations through the Galactic 
Plane are affected by source overlap along the line of sight and 
the optical and UV radiation emitted by young stars is absorbed 
or scattered by the Interstellar Medium (ISM) through the ex- 
tinction process. 

Many SFR values have been measured during the years 
using differe nt datasets and techn i ques, such as infrared (IR) 
photometry llRobitaille & Whitnevi boifi , which measures the 
light from yo ung stars reemitted by the ISM in the IR, free- 
free emission jMurrav & Rahmanl l2010ll . which measures the 
amount of photons required to produce the obs erved ioniza- 
tion of the HII regions, or high-mass star counts [Reed, 2005]. 
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The values r ange from 1 M /y r |Robitaille & Whit nevi 1201 Oil 



to lOM p/yr llGtisten & Mezgertfl982ll . Recently IChomiuk et all 



120 12H observed that, when normalizing all these measurements 
to the same Initial Mass Function (IMF), they converge to 
1.9M Q /yr. 

In recent years, a lot of progress has been made in identi- 
fying the position of many star forming re gions |Russeil, 
Benjamin et al., 2005; Russeil et al., 201 1], hence improving our 
knowledge of the structure of the Galaxy. This has been made 
possible mostly by combining the information from line trac- 
ers with IR surveys, like GLIMPSE [" Galactic Legacy I nfrare d 
Mid-Plane Survey Extraordinaire", [Benjamin et al., 2003], 
MIPSGAL ["A 24 and 70 Micron Su rvey of the Inner Galactic 
Disk with MIPS", ICarev et all l2009ll, and Hi-GAL P 'Herschel 



Infrared Galactic Plane Survey", Molinari et alll2010alfbTl . which 



provide a wealth of data in the domain where the dust sur- 
rounding young stellar objects (YSOs) has the emission peak. 
In these bands it is therefore possible to study the quantity of 
young stars through the UV and optical light of the protostar 
which is absorbed, processed and reemitted by dust in the IR 
domain. In particular, Hi-GAL mapp ed the Galactic Pla ne in 
\{\ < 60 ° and lfrl < 1° with t he PACS BPoglitsch et alll2010tl and 
SPIRE I Griffin et all 1201 (ill instruments in parallel mode, in the 
70, 160, 250, 350, and 500 fim bands. Given its spectral range 
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and sky coverage, it provides a unique opportunity to study col- 
lapsing dust clouds and protostars, i.e. the early stages of star 
formation. Moreover, the survey is particularly sensitive to high 
mass stars (OB, M > 8M ) which regulate the ecology of our 
Galaxy as a whole. 

In this paper, we estimate the SFR in the sky regions mapped 
during the Herschel Science Demonstration Phase (SDP), i.e. 
two 2° ■ 2° tiles centered on the Galactic Plane in t = 30° 
and I - 59°. The t = 30° field observes the Sagittarium and 
Perse us arms, whil e t = 59° is centered on an inter-arm re- 
gion llRusseil et all 1201 111 . Due to its location, the € = 30° 
field is expected to be more active, in terms of star formation, 
than the £ = 59° field. A higher number of YSOs and HII re- 
gions arejocated in the I - 30° area, including the W43 com- 
plex IBal lv et all [T9991. a massive star forming region with an 
associated giant HII region. This gives us the opportunity to 
study the star formation processes both in a very active and in 
a quiescent environment. 

Due to the aforementioned difficulties in estimating the SFR 
of the MW, we do not have a reliable far IR star formation 
indicator. Many indicators have been calibrated on extragalac- 
tic star forming complex es and then applied to Galactic obser- 
vations [see for exam ple iKennicutA Il998t ICalzetti et all l2007t 
iKennicutt et alll2007ll . Since one of the main goals of Galactic 
star formation studies is to derive a SFR estimator of the MW, 
hence test the validity of estimators calibrated o n other galaxies , 
we make use both of an extragalactic indicator IlLi et all [20101 
and of a star counting method to estimate the SFR in the two 
Hi-GAL SDP tiles, and then compare the results. This way, we 
both study the star formation activity in the two Hi-GAL fields 
and test the method in order to apply it to a larger sample when 
the whole Hi-GAL survey will be completed. The paper is or- 
ganized as follows: Sec. [2] describes the dataset; Sec.[3]ouflines 
the SED fitting, the sample selection criteria, and the estimate 
of the candidate YSOs physical parameters such as tempera- 
tures and masses; in Sec. [4] we estimate the SFR with the ex- 
tragalactic estimator and counting the young stars. The approxi- 
mation adopted in this procedure and the quantifiable errors are 
discussed in Sec. [5] Conclusions are summarized in Sec. [6] 

2. The Dataset 



The starting point of this st udy is the catalo g of sources 
of the two Hi-GAL SDP fields llEha et all l20Toll . A new ver- 
sion of this catalog has recently been released to the Hi-GAL 
consortium, obtained with updated data cleaning procedures, 
mapmaking, calibration factors, extraction and photometry in 
the PACS and SPIRE bands. For further information about 
the Hi-GAL pipe line from raw data to map production we re- 
fer the reader to iTraficante et alj 11201 111 . The algorithm used 
for the sou rce detection i s CuTE X (Curvature Thresholding 
Extractor) [Molinari et ajl 1201 ill which double-differentiates 
the sky image and studies the variation of the curvature above 
a given threshold. The identified source profiles are then fit- 
ted with a 2D elliptical Gaussian plus an underlying inclined 
planar plateau. This allows us to detect sources in the presence 
of a variable background, as the one in the Galactic Plane, and 
to select not only point-like sources but compact objects. These 
two conditions are crucial for us because we are interested in 
studying YSOs emitting in the PACS and SPIRE bands, which 
are still embed ded. The considered sources have known kine- 
matic distances llRusseil et all 120 nil and are detected in at least 
three contiguous Herschel bands. This last condition is required 



to exclude spurious detections. After these conditions have been 
applied, we have a total of 681 sources in I — 30° and 316 in 
I = 59°. In order to better sample the SEDs and constrain the 
evolu tionary stage, the 24 //m flux from MIPSGAL llCarev et all 
2009] is also measured at the same position of the 70 /mi 
band. For more information about the sour ce extraction and 
dista nce determination we refer the reader to lElia et~al1 [2010] 
and lRusseil et all ||201 111 . 

3. SED fitting 



The physical parameters of the sources with at least three 
positive fluxes (S \) are estimated through a modified black body 
fit. The emission at wavelength A can be modeled as 



S A (eo, T) = eo 



Ba(T) 



(1) 



where T is the source temperature, eo is the emissivity at the 
reference wavelength Aq = 100 /vm, the emissivity spectral index 
is set to 2 and B A (T) is the blackbody at temperature T. The 
emissivity can be rewritten as 

Mk 



eo = 



d 2 



(2) 



where M is the total mass of the source, Icq is the mass opacity at 
wavelength Aq per unit mass and d is the heliocentric distance. 

We perfor m the fit using a Monte Carlo Markov Chain algo- 
rithm rMCMC Lewis & Brid"lell2002ll and estimate the tempera- 
tures and the tot al masses. We do not fi t the SEDs with the grid 
of models from Robitaille et al. |2006], as was done with the 
previous version of this catalog, because almost all the sources 
peak at waveleng t h A > 160 jum which is the threshold identi- 
fied by lElia et all ll2Q 1 Oil as the one above which the Robitaille 
models do not apply anymore. The fit is performed including 
only Herschel PACS and SPIRE bands. The associated error bars 
come mainly from calibration uncertainties and background re- 
moval. Since we are not analyzing point-like objects but embed- 
ded sources we make use of the calibration errors measured on 
extended emission which amount to 20% of the total flux for 
PACS and 15% of the flux for SPIRE. We also include a 10% 
statistical error coming from the background fluctuations. The 
two errors are added in quadrature and associated to the fluxes. 
In the present analysis we chose not to include the 24 //m flux 
in the fit since we are interested in the envelope emission. The 
24 fj.m flux samples the internal source and this would require 
a two-component model, one for the envelope and one for the 
central star. We do not probe sources with T < 7 K in our anal- 
ysis because, in that temperature regime, the assumption of an 
optically thin regime implicit in the modified black body fitting 
breaks down. We would need a full radiative transfer model to 
account for optical depth effects and this goes beyond the scope 
of the present paper. The average temperature in the two fields 
are < r, 30 >= 17.5 K and < T C59 >= 14.3 K. The I = 59° 
tile is colder on average, consistent with the fact that, being an 
inter-arm field, star formation is not very active. 

In order to estimate envelope masses, we put ourselves in 
the Rayleigh- Jeans (RJ) regime which is optically thin. Since the 
500 /vm band has a low resolution and is affected by source mul- 
tiplicity issues, the 350 /vm band is then the best compromise, so 
we estimate masses according to the formula 



M : 



'350 



&350 ^350(71 



(3) 
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where S 350 is the flux at 350 ji/m, d is the kinematic distance, 
£350 = 0.0 7 cm 2 g~ l is the m ass opacity coefficient at 350 yum per 
unit mass IPreibisc h et al.Lfl993ll . T is the temperature estimated 
by means of Eq. [TJ The measured masses range between 1 and 
10 5 M in € = 30° and between 0.1 and 10 4 M in I = 59°, with 
median values of 454 M and 48 M in i = 30° and I = 59°, 
respectively. 

3. 1 . Selection criteria 

Once we have fluxes, temperatures and masses of the whole sam- 
ple of sources, we perform a further selection to ensure we col- 
lect only YSOs for the SFR estimate. First, we remove sources 
with a negligible emission in the PACS 70 yum band. Second, we 
identify and remove AGB stars using color-color criteria. 

The presence of a 70 yum counterpart has been demon- 
strated to_be_correMedjAdth the internal luminosity of a pro- 
tostar [Dunham et al., 2008] and, moreover, we do not expect 
our survey to detect a 70 fim flux of diffuse ISM just heated by 
the interstellar radiation field. For those reasons, we consider as 
starless all sources without a detection at 70 yi/m. 

One of the most challenging steps in cleaning our catalog is 
to make sure not to include AGB sources in our sample. Around 
30% -50% of the m i d-infra red sources are estimated to be AGBs 
by iRobitaille etall ll2008h using the GLIMPSE survey. AGB 
sources also emit in the IR and can be disentangled from YSOs 
by applying c olor-color criteria. We ma ke use of the criteria de- 
veloped from Mar tinavarro et al.l ll2012ll in the Herschel bands. 
According to these authors, AGBs are located in the color-color 
diagram where the following is verified: 



1=30 



1=59 



0.2 < 



<0.9, 1.2 < 



< 3.3 



(4) 



The percentage of AGBs among sources detected at 70 yum 
is 47% in I = 30° and 52% in t = 59°, which is consistent 
with what was found in the GLIMPSE survey. All the sources 
satisfying these color criteria are removed from the sample. 

In order to identify starless objects which are gravitation- 
all y bound, i.e. prestell ar sources, we apply the method used 
bv lGiannini et al.l Il2012ll on the Vela-C molecular cloud. Instead 
of the virial mass, they make use of the Bonnor-Ebert mass, 
Mbe ~ 2-4-Rbeg 2 /G, where Rbe is the Bonnor-Ebert radius 
which can be approximated by the actual radius of the source, 
Rbe - 6/2 where 6 is the diameter. G is the gravitational con- 
stant and a = -y/^r/yu. Here, Icb is the Boltzmann constant, T is 
the source temperature and fi is the mean molecular weight. The 
source diameters (6) are estimated from the 250 yum map by ap- 
proximating the ellipse provided by CuTEX as a circle with the 
same area and deconvolving it by the instrumental beam. The 
choice of the 250 yum has been done because all sources are de- 
tected in this band. The size distribution is shown in Fig. [1] In 
both tiles the majority of sources h ave 6 > 0.1 pc meaning that 
they are essentially clumps (see e.g. Kauff mann et al.l ll2010ll and 
references therein), as is also confirmed by their mass distribu- 
tion. 

Following the criteria in iGiannini et al.l Il2012ll we identify 
sources with M/Mbe > 0.5 as gravitationally bound. 

The temperature distribution of the prestellar and protostellar 
sources is reported in Fig. [2] This figure clearly shows two pop- 
ulations of objects, a cold one (red histogram), mostly prestel- 
lar, which peaks around 1 1 K and a warm one (blue histogram), 
mostly protostellar, which peaks around 20 K. The protostellar 
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Fig. 1. Diameter distribution of all sources detected in three ad- 
jacent Herschel bands, estimated from the 250 yum map. 
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Fig. 2. Temperature distribution of the considered YSO candi- 
dates and prestellar cores estimated by means of a modified 
black body fit from 70 yum to 500 yi/m. 



sources constitute the final sample used in the following analy- 
sis. 

Tw o recent papers llBattersbv et all 1201 lb iParadis et all 
l2010tl have measured the temperature of these fields on a pixel- 
by-pixel basis, smoothing the signal of all wavelengths to 35.6", 
i.e. the resolution of the Herschel 500 yum band. Their results 
don't show the same bimodal distribution as ours because, in our 
case, we keep the original resolution of each band and, therefore, 
we identify and collect compacts objects and we are sensitive to 
the temperature variations on these smaller scales. 

After applying those criteria, the number of identified candi- 
date YSOs is 235 in I = 30° and 50 in t = 59°. TableQ]gives the 
coordinates, fluxes and distances of the 285 selected sources. 

The distance distribution of the selected sources is shown 
in Fig. [3] In t — 30° (left panel), most of the sources cluster, 
as expected, in the Scutum spiral arm which crosses the line of 
sight twice (at d ~ 5.5 kpc (near) and d ~ 9 kpc (far)) and 
in the Sagittarius arm (at d ~ 3.5 kpc (near) and d ~ 11 kpc 
(far)). The peak centered in d ~ 13 kpc corresponds to sources 
located in the Perseus arm. In t — 59° (right panel) the line of 
sight is tangent to the Sagittarius arm (d ~ 3 - 7 kpc) where 
most of the sources are located. For further information about 
the Galactic distribution of Hi-GAL sources in t he SDP fi elds 
and fo r distance derivation we refer the reader to Russe il et al.l 
112011 . 
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Table 1. List of candidate YSOs in the SDP fields, their MIPS and Hi-GAL photometry and distances 



Source Name 


Glon 


Glat 


F(24yum) 


F(70/im) 


F(160/im) 


F(250/iin) 


F(350ycan) 


F(500/im) 


d 




deg 


deg 


Jy 


Jy 


Jy 


Jy 


Jy 


Jy 


kpc 


G029.0540+0.5062 


29.0540 


0.5062 


0.03 


1.78 


11.68 


11.82 


6.73 


4.27 


8.90 


G029.4191+0.6303 


29.4191 


0.6303 


0.31 


7.36 


7.89 


4.56 


2.96 


1.96 


13.50 


G029. 1367+0.4649 


29.1367 


0.4649 


5.01 


7.77 


14.59 


9.78 


7.41 


2.45 


13.10 


G029. 1612+0.4743 


29.1612 


0.4743 


2.55 


4.33 


18.46 


19.06 


13.34 


6.51 


13.10 


G029. 1581+0.4637 


29.1581 


0.4637 


0.24 


25.74 


32.29 


18.80 


13.96 




13.10 


G029. 1383+0.4195 


29.1383 


0.4195 


0.44 


5.80 


8.15 


3.53 


1.63 


0.97 


14.20 


G029. 1538+0.4161 


29.1538 


0.4161 


0.32 


11.27 


12.53 


4.09 


4.07 




13.30 


G029.8533+0.5491 


29.8533 


0.5491 


0.08 


7.64 


11.83 


6.44 


3.91 


2.64 


4.50 


G029.0238+0.0842 


29.0238 


0.0842 


8.14 


2.30 


18.14 


18.78 


10.33 


8.40 


9.90 


G029.8853+0.4918 


29.8853 


0.4918 


0.05 


1.72 


5.10 


5.71 


5.00 


2.25 


13.40 


G030.4278+0.7694 


30.4278 


0.7694 


0.17 


4.20 


5.54 


2.86 


1.36 


0.50 


11.50 


G029. 1172+0.0889 


29.1172 


0.0889 


44.05 


69.27 


72.11 


41.44 


30.04 


16.74 


9.30 



Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available at this link 
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Fig. 3. Distance distribution of the selected candidate YSOs in 
SDP fields. 



The color-color plots of the final sample are shown in Fig. [4] 
Each axis of these figures report the magnitude difference be- 
tween two bands according to the formula [Ai — Aq] = mi -mo = 
-2.5 log(|^j) where S (A) is the flux at the band A and m is its 
magnitude. The left panel shows a subsample of the right panel 
because some of the protostellar sources are detected at 70 fim 
but not at 24 pm. 

4. SFR estimate 

The SFR in the SDP fields is estimated using two independent 
techniques, one based on an IR extragalactic estimator and the 
other on source counting. Results for the tiles as a whole are 
reported in Tab. [2] 

4.1. Source counts 

Starting from fluxes and distances we estimate the bolometric 
luminosities (Lb„i) of our sources according to the formula 

Lboi = And 2 J ™ S A dA (5) 



Fig. 4. Color-color diagrams for all the protostellar sources in 
the t - 30° (black squares) and t = 59° (red diamonds) fields. 
The color distribution is consistent with embedded sources 
which can be fitted only with a modified black bod y and not 
with an embedded ZAMS model llMolinari et al.Ll2008ll . 



Table 2. SFR estimates 



Method 


Field 




SFRrfZe 

(M /yr) 


SFRmw 

(M G /yr) 


Star counts 


I = 30 
^ = 59 


(9.5 
(1.6 


±4.3) • 10~ 4 
± 0.7) • 10- 4 




IR estimator 


t = 30 


(2.4 


± 0.4) • 10~ 4 


0.71 ±0.13 


^ = 59 


(2.6 


± 1.1) • 10~ 6 


0.10 ±0.04 



where A m j n = 24 yum, A max = 500 yum, d is the heliocentric dis- 
tance and S a is the source flux at the wavelength A. The prestel- 
lar and protostellar source distributions of the L/,,,/ in the Galaxy 
towards the two considered tiles is shown in Fig. [5] 

We then build the M env - h\, a i diagram where M env is the 
envelope mass. This kind of diagrams provides information on 
the evolutionary stage of the source and a prescription of when 
it will join the Zero A ge Main Sequence (Z AMS hereafter) and 
with which final mass IMolin ari et all 120081. The source evolves 
along defined tracks according to the chosen evolutionary model. 
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Fig. 5. Galactic distribution of the bolometric luminosity of 
prestellar (red squares) and protostellar (blue circles) sources in 
the SDP fields. 



evolutionary path, from the beginning of the accelerating accre- 
tion phase to the end of the envelope clean-up phase, i.e. when 
the envelope has disappeared and the object is visible in the op- 
tical bands. They are t f = (2.1, 2.7, 3.5, 4.8, 4.5) • 10 6 yr for 
sources with initial envelope masses in the bins (8 0, 140, 350, 
700, 2000) M Q , respectively llMolinari et al.Ll2008H . Due to Hi- 
GAL sensitivity limits, we are not confident in the detection of 
objects with M em , < 10 M G . These sources are then excluded 
from further analysis and the missing contribution to the SFR, 
coming from the low-mass regime, is estimated in Sec. [5] 

The SFR from candidate high-mass YSOs in each tile is 
obtained by summing up the final masses and dividing them by 
the associated formation time: 



We a dopt the model of c ollapse in turbulence-supported dust 
cores llMcKee et all 1200311 . According to this model, the evolu- 
tion of high-mass YSOs has two main phases. In a first phase, as 
soon as the initial cloud starts collapsing, the luminosity of the 
embedded star increases in an accelerated fashion and the ini- 
tial mass envelope decreases because part of the material feeds 
the central core and part is expelled through molecular outflows. 
The path followed by the protostar in the M env — Lboi diagram 
is almost vertical. At the end of this phase, the star reaches the 
ZAMS, or is very close to it, and is surrounded by an HII re- 
gion. The average time of th is phase is 2.7- 10 s y r for an initial 
envelope mass of 13.5 M Q IMolin ariet ail 120081. In the second 
phase, after the end of the accelerating accretion, the remain- 
ing envelope mass is partially drained by other objects forming 
in the same clump and part of it is expelled through molecu- 
lar outflows. An accretion disk is also present, so part of the 
initial cloud still keeps feeding the central star. The envelope 
clean-up phase ends when the object is visible in the optical 
band and the initial cloud has completely disappeared. Since 
the luminosity in this phase remains almost constant, this results 
in an horizontal track in the diagram. The total average evolu- 
tionary time for the entire path is 3. 5 TO 6 yr for an initial enve- 
lope mass of 13.5 M . For more detailed i nformation about the 
evolutionary mode l we r efer the reader to McK ee et alJ Il2003ll 
and lMolinari et alj Il2008ll . 

Since the sources we are considering are still deeply embed- 
ded, we can assume that M corresponds to the M em , derived in 
Sec. [3] The M em , - Lt, i diagrams for the two tiles are shown 
in Fig. [6] The solid black line is the best log-lo g fit of the high- 
mass counterpart of the low-mass Class I regime [ Mol inari et all 
l2008ll . while the dashed black line is the high-mass counterpart 
of the best log-log fit of the low-mass Class regime. In the high- 
mass domain the subdivision of young stars in classes is a con- 
ventional extrapolation from low-mass stars. Since high-mass 
clumps contain unresolved structures which might fragment into 
more than one high-mass star, and continue to accrete even after 
entering the main sequence, one has to be careful about this sub- 
division. The majority of the sources in our sample are below this 
line, meaning that these objects are still in an early evolutionary 
phase as already argued from the color-color plots in Fig. [4] The 
average uncertainties on luminosities and masses in I — 30° are 

< cr L /L >= 0.2 and < cr M /M >= 0.4, while in in ( = 59° are 

< cti/L >= 0.4 and < cr M /M >= 0.7. These values have been 
estimated through the error propagation of fluxes and distances. 

In order t o estimate th e SFR from the evolutionary tracks 
provided by McKee et al. |2003], we bin these diagrams into 
two-dimensional histograms. Each bin has an associated forma- 
tion time which depends on the source stage and mass. Since 
high-mass objects keep accreting even after they join the ZAMS, 
the formation timescales used in our analysis cover the entire 



Nclasi AW 

SFR = Z Z n M(iJWzAMs(j)/tf(i) M Q /yr (6) 

i=l 7=1 

where Nciass - 3 (counterparts of low-mass Class 0, I, II), 
Nmoss - 5 are the number of classes and of initial masses, 
respectively. riM(i,j) is the number of sources of Class i with 
envelope mass j and MzamsU) is the expected final mass of 
sources of initial envelope mass j when they reach the ZAMS. 
The SFR obtained with this method is (9.5 ± 4.3) • 10~ 4 M Q /yr 
in i = 30° and (1.6 + 0.7)- 10~ 4 M /yr in t = 59°. We do 
not extrapolate these estimates to the whole MW because they 
have been calculated through a counting procedure, so they 
are very local values. Since YSOs generally join the ZAMS 
when the accelerating accretion phase has been completed, 
we report also the SFR estimated at the end of this phase, so 
with shorter evolutionary timescales, in order to provide an up- 
per limit to the estimates. The average evolutionary timescales 
of the main accretion phase for the mass bins mentioned be- 
fore are tf = (4.5, 3.7, 2.7, 2.1, 1.5) • 10 5 yr and the correspon- 
dent SFR would be (1.6 + 0.7) • 10~ 2 M /yr in i = 30° and 
(7.9 ± 3.6) ■ 10~ 4 M /yr in t = 59°. 

4.2. Monochromatic estimator at 70 \im 

Another way to study the SFR of a sample of sources with known 
infrared luminosities is to make use of an extragalactic estima- 
tor. Dust surrounding the forming stars absorbs UV radiation 
and re-radiates it in the IR. Therefore, optical/UV based indica- 
tors are not reliable in the MW because of the high extinction 
generated by interstellar dust along the line of sight. On the con- 
trary, IR emission is a reliable tracer of th e photons produced by 
the forming source. lLawton et al.1 Il2010ll finds that, among all 
the IR bands, the 70 fim band is the most reliable indicator of 
star formation because it is abl e to account for IR emission from 
HII r egions [see for example iTibbs et ail 1201 2t iFaimali et all 
1201211 . 

We the n consid e r the 70 jum monochromatic estimator de- 
veloped by iLi et al. I i20loll . They measure the 70 /jm luminosi- 
ties of extragalactic star forming regions and calibrate them on 
their previous SFR estima tor based on 24 fim and H(a) luminosi- 
ties llCalzetti et all 1200711 . The conversion factor between total 
luminosity at 70 /zm and the SFR is 1.067 ■ 10 43 erg/s. The SFR 
obtained with this method is (2.4 ± 0.4) • 10~ 4 M G /yr in I = 30° 
and (2.6 ± 1.1)- 10~ 6 M /yr in ^ = 59°. 

Since the estimator is built averaging extragalactic star form- 
ing complexes, we extrapolate these values to the whole Galaxy 
modeling it as a disk in which sources are uniformly distributed 
within the volume. The observed fields are then slices of this vol- 
ume, with the vertex centered on the Sun and an aperture angle 
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Fig. 6. M em ,-Lt, i diagram of the sources in the £ = 30° (left panel) and £ = 59° (right panel) fields for the prestellar and protostellar 
cores shown in Fig. [2] The solid black line and the dashed bla ck line are the best lo g-log fit of the high-mass counterparts of the 
Class I and Class sources in the low-mass regime, as found in lMolinari et al] ll2008ll . respectively. 



of 2°, i.e. the dimension of the tiles. The SFR value of the single 
tiles is then extrapolated to the whole MW through the equation: 



1=30 



1=59 



SFRmw — SFR t 



where 



Vmw 

~v7 



Vm 



= nR l h h 



(7) 



(8) 



dfa 
V, = -L-t, 

In the previous equations, V t and Vmw are the volumes of 
the tile and of the whole MW, respectively. Iimw = 0.3 kpc is 
the thickness of the Galactic disk and R = 15 kpc is its radius, 
d t is the average heliocentric distance of the candidate YSOs in 
the tile (7.6 kpc in t = 30° and 3.3 kpc in i = 59°), a = 2° 
is the aperture angle and Zt i s the distance of the sources from 
the Galactic plane. Following Pala dini et al.l ll2004ll . we estimate 
Zt by means of a Gaussian fit of the z = dsin(b) distribution 
of the sources, where b is their Galactic latitude. The FWHM 
of the fitted Gaussian curve is the z, value we assume for all 
the protostars in the same tile. The fits are shown in Fig. Q 
We obtain zm — 71.9 pc and zm - 27.6 pc. These values are 
in agreeme nt with the di stribution of 456 Galactic HII regions 
reported in Paladini et al.1 Il2004ll . 

We obtain a SFR W = 0.71 +O.13M /yr in I = 30° and 
SFR M ,y = 0.10 + 0.04 M G /yr in I = 59°. As expected, the 
extrapolation from t = 30° gives more reasonable results 
since it is an area with active star forming regions and, there- 
fore, it has characteristics similar to the star formation domi- 
nated galaxies used to calibrate the estimator. In fact, the value 
of SFRmw from I = 30° is in good agreement with previ- 
ous resul ts in the MW obtaine d using different estimators and 
datasets IChomiuk et all 1201 2l| and it falls in the range mea- 
sured by Rob itaille & Whitney] ll2010ll using the Spitzer/IRAC 
GLIMPSE survey (0.68 - 1.45 Af e /yr). In the previous extrap- 
olation of the local SFR to the whole Galaxy we set the overall 
tile distance d, (Eq. [8]) to the average distance of the sources in 
the considered tile. When we set d, to the distance of the further 
source in the tile (d t = 14.4 kpc in I — 30° and d t = 7 kpc in 
t = 59°), i.e. the distance in which we assume our sample to be 





-100 -50 50 100 

z (pc) 



Fig. 7. z distribution of the Hi-GAL protostars in the I = 30° and 
i = 59° fields. The z value adopted to estimate the tile volumes 
is the FWHM of the Gaussian fit. The second peak at z ~ 60 pc 
in £ — 59° is d ue to a concentrati on of YSOs in the Vulpecula 
OB association [Billot et alll2010ll . 



complete, we obtain SFR MW = 0.20 + 0.04 M /yr in i = 30° and 
SFRmw = 0.02 ± 0.01 M /yr in £ = 59°. 



5. Catalog completeness 

The catalog completeness provides us with the information 
about how many sources we can miss with fluxes fainter than 
a given threshold. This translates into how much mass and lu- 
minosity we are missing when calculating the SFR. Since the 
masses have been estimated from the 350 yum flux (Eq. 0), we 
study the distribution of the fluxes in this band to identify the 
lower detection threshold. In Fig. [8] we show a zoom of the 
70 fim (top line), 160 fim (central line) and 350 fim flux dis- 
tribution (bottom line) of the candidate YSOs in the two tiles. 
We can see that below a given flux value less and less sources 
are detected and that this threshold is lower in I — 59° than 
in I — 30°. The threshold value depends mostly on the instru- 
mental sensitivity, on the detection algorithm and on the ob- 
served field. The more the field is populated, as in the i = 30° 
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case, with bright sources and variable background, the more dif- 
ficult it is for the detection algorithm to detect faint objects. We 
chose as a flux sensitivity threshold at 350 fim the values right 
before the first peak of sources. These limit fluxes are 1.5 Jy 
and 0.7 Jy for ( = 30° and t = 59°, respectively. Thes e val- 
ues are very similar to the ones estimated by iMolinari et alj 
ll2010all through synthetic source experiments in the old ver- 
sion of the catalog. In order to translate those values into the 
amount of missing mass, we estimate the mass corresponding 
to those limit fluxes using Eq. [3] where we assign, as distance, 
the median distance detected in the field and, as temperature, 
the average temperature of protostars in that field. The conser- 
vative choice of the median distance associated to the limit flux 
allows us to include even the faintest objects in the estimate of 
missing mass. The envelope masses corresponding to the limit 
flux then turn out to be M"™, = 66 M Q and M 1 ™, _ 6.5 M in 
I = 30 ° a nd I _ 59°, respective ly. According to Molinari et alj 
|2008| and lSaraceno et al] fl996], this translates into a final stel- 
lar mass of M<™ = 6.5 M Q and M l ™ = 0.8 M . The nu mber of 
source s expected with M < Mum is estimated with the iKroupal 
11200 111 Initial Mass Function (IMF). The IMF is normalized to 
the counts in the bin of lowest envelope masses (80 M < 
AMq < 140 M ) where our detection is complete, according 
to the formula 

iV(M.),W— |-] (9) 

where N(M„,j„) is the number of stars in the bin AM m ,„ = 
[0.1, 0.5] M , N(M ) is the number of stars in the bin AM and 
Mo is the limit mass of the final star M hm in each tile, y = 2.3 
for M > 0.5 M and y = 1.3 for 0.1 M < M < 0.5 M . We 
then calculate the number of missing stars until M„„„ =0.1 M , 
which we assume to be the minimum value for a YSO mass. By 
means of this procedure we estimate that, considering 0. 1 M as 
lower limit for a final star mass, we are not detecting ~ 55% of 
the total number of sources in i — 30° and ~ 3% in t — 59°. 
The average timescale for a low-mass YSO population i s t ~ 2 
Myr (see fo r example lEvans et al.1 ll2009ll : ICovev et all ll2010ll : 
lLadaetal.1 Il2010ln . By using again the star counts method, we 
can then calculate the amount of missing SFR as SFRM<M(iim) = 
M min •A^(M„„„)/r. We obtain SFR M <M(/ lm ) = 1 .0- 10~ 5 M /yr and 
SFR M<M(/ , m) = 3.9 • 10~ 8 M /yr in i = 30° and t = 59°, respec- 
tively. Those values are negligible with respect to the estimates 
of the SFR obtained with both methods. 

6. Summary and conclusions 

We estimate the SFR in the two best studied Hi-GAL fields mak- 
ing use of the most updated version of the source catalog. Due to 
the spectral range, sensitivity and sky coverage of our dataset and 
because of the extraction algorithm, we are particularly sensitive 
to the very early stages of high-mass star formation. Moreover, 
the two Hi-GAL fields in this study cover both Galactic arms 
and inter-arm regions and this gives us the chance to study the 
SFR both in a very active and in a more quiescent field. Thanks 
to color-color criteria and to the knowledge of heliocentric dis- 
tances, we are able to remove AGBs stars, keeping only the pro- 
tostellar sources. 

In this work we make use both of an extragalactic monochro- 
matic estimator at 70 fim and of a source counting procedure 
and compare the results. The SFR estimates for the two tiles 
are (9.5 + 4.3) ■ 10~ 4 M /yr and (1.6 + 0.7) • 10~ 5 M /yr with 
the star counts method in t — 30° and { - 59°, respectively. 



1=30 1=59 
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Fig. 8. Distribution of the candidate YSOs 70 fim, 160 fim and 
350 fim fluxes in the SDP fields. These bands are particularly 
important because the 70 fim fluxes are used as indicator of star 
formation activity, they 160 fim is located in the peak of the SED 
and the 350 fim fluxes determine the masses through Eq. [3] The 
vertical red-dashed lines in the bottom panels indicate the sen- 
sitivity limit used to calculate the completeness of our sample. 
They correspond to 1.5 and 0.7 Jy in I = 30° and € = 59°, 
respectively. 



Results with the IR estimator are (2.4 + 0.4) • 10~ 4 M /yr and 
(2.6 + 1.1) • 10~ 6 M /yr in t = 30° and t = 59°, respectively. 
Those values are in good agreement mi- 30°, which is an 
active field. 

We also estimated that our catalog is complete at ~ 45% in 
I = 30° and ~ 97% in I = 59°. Because of the instrument char- 
acteristics, ml- 30° we miss a large population of low-mass 
stars while the i = 59° field is closer and less populated, so the 
sampling is more complete. The contribution to the Galactic star 
forming activity coming from non detected low-mass stars is es- 
timated to be at least more than one order of magnitude lower 
than the SFR from high-mass objects obtained using both meth- 
ods. 

The extragalactic indicator is, by construction, an average 
value of extragalactic star forming regions, so it provides SFR 
estimates which are meant to be calculated on whole galaxies. 
Therefore, it makes sense to extrapolate the local results of the 
tiles to the whole MW. When applying this method to the i - 30° 
tile, which has a similar star formation activity as the galaxies 
on which the indicator has been calibrated, we get results in 
agreement with previous studies (S FRmw ~ 0.71 +0.13 M /yr) 
and the two methods provide consistent results. This means that, 
when applied to areas with characteristics similar to the extra- 
galatic star forming regions used to calibrate it, the IR estimator 
gives reliable results even in the MW. On the contrary, the SFR 
from the source counting method is not meant to be extrapolated 
because is based on the local population of YSOs, so it would 
provide very different results depending on the observed area. 

Both methods are dominated by model errors which are 
not quantifiable. The evolutionary star model, the IR indicator 
and the Galactic extrapolation are subject to assumptions which 
might be negligible on averages of large numbers but might be 
significant when a limited sample is available. Therefore it will 
be important to expand the analysis described in this work to 
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the whole Hi-GAL survey, which will cover the entire Galactic 
Plane and, with it, the vast majority of the star formation com- 
plexes. The use of a significantly larger sample of YSOs dis- 
tributed across the Galaxy will also allow us to derive an effec- 
tive Galactic SFR indicator. 
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